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Changes  of  the  Sq  quiet  daily  variations  of  the  geomagnetic 
field  as  a  function  of  time  of  year  and  latitude 
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Abstract 

'An  analytic  representation  of  the  quiet  daily  variation  of  the 
geomagnetic  field,  Sq,  was  obtained  for  the  North  American  region  with  a  mean 
error  of  less  than  10Z  of  the  daily  variation  range-  The  month  by  month  24-, 
12-,  8-,  and  6-hr  Fourier  components  of  the  field  in  H,  D,  and  Z  were  found  to 
have  significant  annual  and  semiannual  changes  that,  in  general,  agreed  with 
the  expectations  for  ionospheric  current  sources  for  the  changes.  Tables  and 
computer  programs  are  provided  from  which  the  quiet  field  representation  may 
be  constructed  for  any  selected  North  American  location.  At  the  equator  the 

Sq(H)  variation  was  over  125  gammas  and  showed  significant  equinoctial 

-J,->  *1  .J 

maxima.  At  30°  to  40°  the  lowest  values  of  Sq(H)  (less  than  20  gammas)  were 
representative  of  the  region  of  the  ionospheric  current  vortex.  Mid  latitude 
Sq(D)  values  were  about  50  gammas  in  range.  Great  enhancements  of  the  daily 
field,  particularly  in  Sq(H)  and  Sq(Z)  occurred  at  high  latitude  with  a  strong 
summertime  maxmura  range  of  about  90  gammas  in  Sq(H)  at  about  75°  lat.  Five- 
degree  latitude  representations  of  the  quiet  daily  field  change  illustrate  the 
variations  occurlng  during  the  day  and  in  different  months  of  the  year. 
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I-  Introduction 


On  days  that  may  be  described  as  "quiet"  for  the  solar  terrestrial 
environment  of  fields  and  particles  the  magnetograras  of  the  geomagnetic  field 
variations  at  the  earth's  surface  show  regular  daily  changes  of  the  field  that 
are  a  few  tenths  to  several  percent  of  the  strength  of  the  earth's  main 
field.  Figure  1  illustrates  such  a  Sq  variation  in  one  component  of  the 


Figure  1.  Near  here 

field  at  a  mid-latitude  location.  The  figure  also  shows  the  concurrent  zenith 
angle  x  of  the  sun  and  (Cosx)  called  the  Chapman  factor,  that  is 

approximately  proportional  to  the  ionospheric  Ionization  near  100  km  for  the 
day.  The  solar  control  of  the  quiet  daily  variation  Is  immediately  evident  in 
the  enhancement,  onset,  and  cut-off  times  of  the  0  component  of  the  field. 

Any  spectral  analysis  of  the  day's  variation  of  the  magnetic  field  at  a 
station  on  the  earth's  surface  shows  clearly  defined  peaks.  Figures  2 a  and 


Figure  2.  Near  here 

2b  Illustrate  a  typical  spectral  composition  with  peaks  of  constantly 

decreasing  amplitude  at  24,  12,  8,  and  6  hr  and  the  absence  of  outstanding 

features  at  shorter  periods.  These  spectral  peaks  that  arise  because  of  the 

dependence  of  Sq  upon  solar  time  provide  a  convenient  way  for  decomposing  the 

* 

geomagnetic  dally  variation  into  analytic  parts  for  detailed  study. 

The  Sq  pattern  at  a  station  undergoes  month  to  month  changes*  The  12- 
month  display  In  Ptgure  1  Illustrates  how  the  maxima  and  minima  in  field 
excursions  can  change  In  time  and  amplitude  over  the  year.  Both  annual  and 
semiannual  variations  are  evident  here.  Such  variations  occur  In  many 
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geophysical  phenomena  of  ionsopheric  and  magnetospherlc  origin. 

The  Sq  pattern  changes  with  latitude.  Figure  3  shows  how  different  the  H 


Figure  3-  Near  here 

component  of  field  at  Huancayo  is  from  the  same  component  at  Fredericksburg  or 
Toolangi  and  how  the  annual  variation  in  time  of  minimum  differs  between  these 
last  two  stations  that  are  in  different  hemispheres.  The  magnetograms  in  the 
21-June  column  of  Figure  4  indicates  how  the  quiet-day  variation  will  appear 


Figure  4.  Near  here 

at  a  number  of  locations  from  the  equator  to  the  pole- 

When  the  magnetic  field  is  active  (the  15-June  column  of  Figure  4)  the 
station  magnetograms  may  bear  little  resemblance  to  those  on  the  well-behaved 
(i.e.,  consistent,  slowly  changing)  quiet  days.  The  very  best  pictures  of 
earth's  quiet  field  behavior  are  obtained  by  a  selection  of  data  from  those 
years  that  have  experienced  the  lowest  number  of  solar-terrestrial 
disturbances.  Such  occasions  happen  approximately  every  11  years  as  a  result 
of  the  solar  activity  cycles. 

If,  from  each  day's  geomagnetic  record,  the  average  of  the  field 
strengths  for  the  quietest  days  of  the  corresponding  month  are  subtracted,  the 
residuals,  when  averaged  in  lunar  time,  will  show  a  lunar  semidiurnal 
variation  in  phase  with  the  atmospheric  tides.  (Lunar  time  is  computed  for  a 
24-hr  duration  between  times  of  largest  lunar  elevation  angles,  called  lunar 
noon).  Figure  5  shows  the  lunar  field  variations  in  D  component  at  Honolulu 
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Figure  5.  Near  here 

together  with  the  Fourier  component  sine  waves  obLaLned  from  the  spectral 
analysis  of  this  variation.  The  semidiurnal  variation  Is  the  most  significant 
component.  Figure  6  illustrates  how  much  this  lunar  component  contributes  to 


Figure  6.  Near  here 

the  observed  field  at  a  midlatitude  station  a  particular  day. 

The  present  study  of  Sq  will  be  limited  to  a  description  of  the  4 
principal  spectral  components  (24,  12,  8,  and  6  hr)  as  a  function  of  latitude 
with  the  annual  and  semiannual  changes  included.  Only  the  North  American 
quadrant  of  the  earth  will  be  considered  for  this  representative  pilot 
program.  The  quietness  condition  was  best  met  in  1965,  a  year  of  extremely 
low  geomagnetic  activity  from  solar-terrestrial  sources.  The  daily,  3- 
component  (H,  D,  and  Z)  geomagnetic  recordings  will  be  separated  by  months 
from  which  days  of  particularly  low  activity  will  be  selected  as 
representative  of  Sq.  Lunar  changes  will  be  removed.  Using  standard  spectral 
analysis  techniques  on  the  data  from  15  observatories  the  Sq  field  variation 
will  be  represented  analytically  as  a  function  of  latitude.  A  short  computer 
program  will  be  designed  to  reconstruct  the  daily  field  change  at  any 
location.  An  error  analysis  will  determine  the  adequacy  of  such  a  program  to 
represent  the  field  values  at  the  observatories.  At  the  present  stage  of  the 
study  no  attempt  will  be  made  to  separate  the  Induced  Internal  sources  from 
the  external  sources  of  the  measured  fields. 
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II •  Background 

Before  I  describe  the  analysis  and  results  let  me  review,  briefly,  what 
is  known  of  the  source  of  the  Sq  variations.  The  best  interpretation  of 
Sq-field  changes  supports  a  model  that  pictures  the  dynamo  effect  of  diurnal 
and  semidiurnal  tidal  winds  in  the  lower  ionosphere  generated,  in  situ,  by  the 
solar  heating  (Richmond,  1979)  •  The  "dynamo  effect"  is  a  term  used  to 
describe  the  physical  process  like  that  occurring  in  a  hydroelectrical  plant 
to  convert  energy  of  mechanical  rotation  into  eletricity.  The  thermally 
driven  atmospheric  winds  transport  a  conducting  ionosphere  through  the  earth's 
main  field  and  thereby  a  current  is  caused  to  flow  in  a  direction  defined  by 
that  field  and  the  driving  wind.  The  ease  with  which  the  current  flows 
depends  upon  the  conductivity  of  the  ionosphere,  a  gas  of  charged  particles 
which,  unlike  most  metal  conductors,  has  a  conductivity  that  depends  upon  the 
density  of  the  various  charges  and  how  they  are  allowed  to  move  in  the 
presence  of  the  earth's  field.  A  daily  change  in  the  electron  and  ion  density 
of  the  atmosphere  is  brought  about  by  the  exposure  to  the  sun  as  the  earth 
rotates.  Figure  1  showed  how  the  solar  zenith  angle  x  and  ionization  factor 
(Cosx)*^  changes  throughout  a  typical  day.  In  Figure  7  the  (Cosx)^^  is 


Figure  7.  Near  here 

displayed  as  a  function  of  latitude  for  the  summer  solstice,  winter  solstice 
and  equinox.  Note  that  the  summer  to  winter  difference  will  give  an  annual 
variation  in  ionization  at  all  locations  except  the  equator  and  a  most  extreme 
annual  variation  at  the  high  latitudes.  There  is  a  large  semiannual  variation 
of  (Cosx)^^  near  the  equator  where  the  equinoctial  values  exceed  the 
solstitial  values  because  of  the  semiannual  overhead  passage  of  the  sun. 
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Figure  8  illustrates  the  seasonal  variation  of  \  and  (Cosx) '  at  a 


Figure  8.  Near  here 

mid-latitude  location.  A  clear  annual  component  in  the  ionization  is  seen 
here.  Also  note  the  asymmetry  between  the  maximum  in  the  Chapman  ionization 
factor  for  the  summer  and  winter  with  respect  to  the  equinoxes;  such  an 
asymmetry  will  necessarily  introduce  some  semiannual  component  into  the  Sq 
analysis.  The  ionospheric  winds  associated  with  the  atmospheric  circulation 
patterns  show  the  annual  and  semiannual  changes  in  solar  heating  and  in 
addition  they  reflect  the  annual  and  semiannual  variations  of  changing  auroral 
activity.  Figure  9  shows  the  mean  meridional  circulation  in  the  earth's 


Figure  9.  Near  here 

thermosphere  for  different  activity  levels-  There  are  not  too  many  direct 
measurements  of  the  wind  systems  that  drive  the  ionospheric  currents;  the 
atmosphere  is  too  dense  there  for  satellite  observations,  the  few  rocket 
measurements  that  have  been  produced  are  relatively  Instantaneous  samples,  and 
the  meteor-trail  studies  are  limited  in  time  and  number.  Global  pictures  of 
the  ionospheric  winds  are  often  Inferred  from  the  Sq  current  Itself. 

Near  the  equator  the  earth's  main  field  is  parallel  to  the  earth's 
surface  at  ionospheric  altitudes.  This  feature,  combined  with  the  more 
Intense  solar  ionization  and  heating,  leads  to  a  unique  equatorial  enhancement 
of  the  Sq  current  system  called  the  "equatorial  eletrojet."  The  influence  of 
this  electrojet  disappears  at  15°  to  20°  from  the  equator. 
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Models  of  Ionospheric  current  patterns  associated  with  Sq  show  foci  near 
noon  at  about  30°  north  and  south  latitude  (Matsushita,  1967).  The  existence 
of  current  foci  means  that  the  ratio  of  the  northward  to  eastward  components 
of  field  changes  rapidly  with  latitude  near  that  location.  Tarpley  (1973) 
showed  a  seasonal  variation  In  the  location  of  these  foci  that  had  a  large 
annual  component  with  conspicuous  harmonics. 

The  semidiurnal  lunar  tidal  oscillations  of  the  atmosphere  drag  the 
Ionosphere  through  the  earth's  field  producing  another  ionospheric  dynamo 
current  that  is  also  dependent  upon  both  the  existence  of  day-side  Ionization 
and  the  field  direction.  Special  analysis  of  this  geomagnetic  lunar 
variation,  L,  (Campbell,  1980)  makes  use  of  the  fact  that  the  lunar  tidal 
"day"  Is  about  50.5  min  longer  than  the  solar  day.  Although  L  Is  about  20 
times  smaller  than  Sq  at  most  mid-latitude  locations  Its  contribution  to  the 
observed  dally  change  Is  appreciable  near  the  equator  and  poles.  A  proper 
analysis  of  the  global  characteristics  of  Sq  should  have  the  L  removed. 

The  earth's  dipole  field  extending  into  space  Is  distorted  by  the 
prevailing  wind  of  field  and  particles  blown  outward  from  the  sun.  Figure  10 


Figure  10.  Near  here 

shows  the  contours  of  the  average  change  from  the  dipole  field  in  the 
magnetospherlc  regions  near  the  earth  at  quiet  times.  Although  we  consider 
the  dipole  field  to  be  locked  to  the  earth  as  it  rotates,  a  surface  location 
will  move  under  the  magnetospherlc  distortion  of  the  type  shown  In  Figure  10 
giving  a  Sq-type  variation  not  of  ionospheric  origin.  As  the  earth's  postlon 
with  respect  to  the  sun  changes  through  the  year  we  may  expect  an  annual  and 
semiannual  change  In  this  magnetospherlc  component  of  the  surface  field 


7 


measurement*  For  example.  Figure  il  shows  the  seasonal  variations  In  the 
midnight  field  levels  observed  at  low  latitudes  which  are  consistent  with  the 
distorted  nightslde  magnetosphere  (Campbell,  1981).  There  Is  still  some 
disagreement  among  geophysicists  regarding  the  field  amplitudes  at  the  earth's 
surface  be  expected  from  the  raagnetospherlc  changes.  It  Is  Important  to 
realize  that  the  sharp  onset  and  cutoff  of  Sq  at  Ionospheric  sunrise  and 
sunset  (c.f..  Figure  1)  cannot  be  ascribed  to  any  magnetospheric 
characteristics. 

At  the  polar  edge  of  the  auroral  zones  the  satellite  evidence  shows  that 
there  are  luminous  auroras  even  during  the  quiet  times  of  solar-terrestrial 
activity.  Such  evidence  means  that  there  exists  an  enhanced  ionization  and 
conductivity  at  this  particle  precipitation  region  and  some  field-aligned 
currents  between  the  magnetosphere  and  the  ionosphere.  The  associated  field 
changes  generally  have  maxima  In  the  midnight  to  dawn  hemisphere.  The  auroral 

zone  ionospheric  currents  connect,  via  the  Ionosphere,  through  the  auroral 

zone,  across  the  polar  cap  and  to  subauroral  latitudes.  The  aurorally  related 
disturbances  display  maxima  at  the  equinoxes.  The  fluctuation  of  auroral 

currents  with  solar  activity  should  give  an  increased  error  In  the 
determination  of  high  latitude  Sq  on  even  relatively  quet  days. 

For  the  polar-cap  observatories,  when  the  average  Sq  quiet-day  field 
variation  is  subtracted  from  the  records,  the  residual  shows  a  variation, 
particularly  evident  in  Z,  that  has  been  associated  with  the  direction  of  the 
solar  wind  interplanetary  magnetic  field  toward  or  away  from  the  sun 
(Campbell,  1976).  Because  this  field  direction  Is  determined  as  the  wind 
leaves  the  sun's  close  environment,  the  earth's  polar  region  response  is 
called  the  "solar  sector"  effect.  The  surface  manifestation  of  this  effect  is 
probably  due  to  the  small  shift  in  position  of  the  polar-cap  quiet  time  Sq 


currents  responding  to  the  magnetospherlc  tail  as  it  interacts  with  the 
parallel  solar  wind  field.  In  deriving  a  mean  monthly  quiet-day  field 
measurement  in  the  polar  cap  a  balance  of  toward  and  away  sector  days  need  be 
included.  Because  this  balance  is  not  always  obtained  the  sector  effect  adds 
some  error  to  the  Sq(Z)  determinations  at  the  polar-cap  observatories. 

Figure  12  is  a  diagram  of  the  complex  array  of  items  that  can  contribute 


Figure  12.  Near  here 

to  the  ionospheric  current  system.  One  objective  of  the  Sq  research  is  to 
determine  the  relative  importance  and  predictability  of  these  contributors. 
In  the  present  study  I  will  try  to  represent  the  quiet  day  variations  in  an 
analytic  form  with  a  limited  number  of  terms  and  to  display  the  expected  quiet 
field  changes  in  North  America  as  a  function  of  latitude  and  local  time.  I 
will  try  to  specify  the  annual  and  semiannual  changes  in  a  way  that  will  allow 
a  comparison  of  the  ionospheric  and  magnetospherlc  sources  of  this  change. 

III.  The  Data  and  Analysis  Technique 

A  large  number  of  the  world's  1965  geomagnetic  observatory  records  have 
been  digitized  at  2.5  min  intervals  and  are  available  on  computer  tapes  at  the 
World  Data  Center  A  for  Solar-Terrestrial  Physics  at  Boulder,  Colorado. 
Because  this  particular  year  was  an  unusually  quiet  one  the  geomagnetic 
records  can  provide  a  remarkably  complete  picture  of  the  quietest  conditions 
of  the  earth's  magnetic  field.  From  these  archives,  data  tapes  were  selected 
that  gave  digital  values  of  the  northward  (H) ,  eastward  declination  (D) ,  and 
vertical-downward  (V)  field  components  for  the  15  observatories  listed  in 
Table  1.  This  sample  of  locations  represented  a  mean  longitude  position  of 
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Table  1.  Near  here 

about  98.5°  West.  The  first  two  stations  of  this  table  are  considered  to  be 
at  equatorial  locations;  the  next  six,  at  middle  latitudes;  the  following 
four,  at  the  auroral  region;  and  the  last  three,  at  the  polar  cap. 

The  station  data  were  arranged  in  geomagnetic  latitude.  Dip  latitude  is 
generally  more  effective  in  the  equatorial  region  because  the  electroject 
phenomena  there  are  controlled  by  field  directions  In  the  local  ionsosphere. 
However,  for  the.  American  hemisphere,  the  differences  between  dip  and 
geomagnetic  equators  Is  negligible  with  respect  to  the  scale  of  analysis  in 
this  present  study.  High  latitude  phenomena  are  better  ordered  in  geomagnetic 
coordinates . 

The  D  components  of  field  had  been  recorded  as  angular  declinations. 
These  values  were  converted  to  eastward  field  strengths  in  gammas  (1  gamma  ■ 
l0-<*  Tesla)  using  the  usual  relationship:  U(gammas)=H0  s ln(D(degrees) ) ,  where 
Hq  is  the  total  horizontal  field  strength  in  gammas. 

The  observatory  data  were  separated  by  month  and  into  directional, 
components.  For  each  month  the  quietest  UT  days  were  selected  by  the 
condition  that  for  all  8  of  the  3-hr  index  intervals,  the  planetary 
geomagnetic  activity  index,  Kp  was  less  than  or  equal  to  2+.  This  is  the 
lowest  value  of  Kp  that  would  provide  a  sufficient-size  sample  of  quiet  days 
in  each  month  of  this  analysis.  For  each  of  the  2.S  min  data  samples  of  those 
days,  the  average  values  were  formed  for  the  H,  D,  and  7  components  separately 
and  called  Sq .  Obvious  erratic  changes,  such  as  that  near  5  UT  on  June  21  at 
Baker  Lake  in  Figure  4  were  removed  from  the  data.  The  appropriate  lunar 
tidal  component  (Campbell,  1980)  was  then  subtracted  and  the  residual  was 
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Table  1.  Observatory  scions 


Station 

„  ,  UT  of 

Code  local  MDT 

Geographic 

degrees 

Geomagnetic 

degrees 

Latitude 

E. 

Longitude 

Latitude 

E. 

Longitude 

lluancayo 

HU 

5.02 

-12.05 

284.67 

-0.64 

354.27 

Tatuoca 

TT 

3.23 

-  1.20 

311.48 

9.51 

21.25 

Paramaribo 

PA 

3.68 

5.82 

304. 78 

16.93 

14.83 

Honolulu 

HO 

10.53 

21.32 

202.00 

21.17 

266.99 

San  Juan 

SJ 

4.41 

18.12 

293.85 

29.57 

3.63 

Tucson 

TU 

7.  39 

32.25 

249.17 

40.48 

312.72 

Fredericksburg 

FR 

5.16 

38.20 

282.63 

49.54 

350.42 

Victoria 

VI 

8.23 

48.52 

236.58 

54.27 

293.51 

Sitka 

SI 

9.02 

57.07 

224.67 

60.09 

275.86 

College 

CO 

9.86 

64.87 

212.17 

64.73 

256.99 

Barrow 

BW 

10.45 

71. 30 

203.25 

68.64 

241.55 

Baker  Lake 

BL 

6.40 

64.33 

263.97 

73.81 

315.98 

Resolute  Bay 

RB 

6.33 

74.  70 

265.10 

83.07 

290.20 

Alert 

AT 

4.17 

82.50 

297.50 

85.87 

167.29 

Thule 

TH 

4.61 

77.48 

290. 83 

88.95 

3.56 

Note:  For 

latitude 

grouping  purposes, 

the 

first  two  stations  of  this 

list 

are  considered 

to  be  at 

the  equatorial  locations;  the  next 

six,  the  middle 

latitudes;  the 

following 

four,  the 

auroral 

regions  ;  and  the 

last  three, 

the 

polar  cap.  UT 

is  given 

in  decimal 

hours . 
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annual  mean  of  the  24-hr  Fourier  component  of  Sq(H)  at  20°  latitude  has  an 
amplitude  of  9.6  gammas  that  occurs  at  10:30  a.m.  local  time. 

Next,  for  the  annual  variation  part  of  Sq ,  the  amplitude  (gammas)  and 
phase  (hr  of  maximum)  are  separated  In  the  tables  to  show  the  annual  change 
(gammas)  of  amplitude  from  the  mean  value  and  the  "month”  of  occurrence  of  the 
maximum  In  this  amplitude  change.  Following  this  Is  the  change  In  hours 
(decimal)  of  the  mean  phase  and  "month"  of  occurrence  of  the  maximum  In  this 
phase  change.  Here,  the  "month"  is  represented  by  a  division  of  the  year  into 
12  equal  parts,  with  January  1  representing  0.  A  repeat  of  this  pattern  of 
presentation  for  the  semiannual  variations  completes  the  six  basic  parts  of 
the  24-hr  component  presentation.  The  tables  then  continue  with  similar 
representations  for  the  12-,  8-,  and  6-hr  field  components. 

In  all,  40  values  are  used  to  represent  the  Sq  quiet  dally  variation  of 
one  directional  component  at  one  latitude.  In  order  that  the  values  for  every 
5°  be  shown,  a  40  x  19  matrix  Is  required  for  each  of  the  three  components. 
The  analytic  representation  of  the  Sq  by  the  coefficients  of  the  three  tables 
provide  a  rather  complete  picture  of  the  behavior  of  Sq  in  North  America,  as  1 
shall  show  shortly.  For  an  example  of  reading  values  directly  for  the  table, 
let  us  ask  how  the  amplitude  of  the  8-hr  component  of  Sq(H)  changes  In  Its 
annual  variation  at  40°  latitude.  From  the  table  we  read  that  there  is  only 
0.2  gammas  change  and  that  this  change  reaches  Its  maximum  In  mid  August 
("month"»7 .6) . 

A  first  picture  of  the  latitude  change  of  the  Sq  amplitudes  can  be 
obtained  from  a  comparison  of  values  along  a  row  in  the  tables.  For  example, 
the  annual  mean  amplitudes  in  H  for  the  24-hr  component  (top  row  of  Table  21 
varies  from  39.0  gammas  at  11:12  a.m.  at  the  equator  to  a  minimum  of  4.1 
gammas  at  5:30  a.m.  at  35°  latitude,  then  to  an  85°  maximum  of  26.8  gammas  at 
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9  p.m. 


The  next  step  in  the  processing  of  the  data  was  t  ■>  produce  a  computer 
program  that  would  plot  the  expected  9q  quiet  datlv  variations  for  anv  desired 
day  In  the  year  and  any  location  In  North  America,  given  the  Tables  2,  1,  and 
4  as  reference  files.  The  Equation  Ml  uses  the  table  *  irtn  of  the  Fourier 
coefficients  to  reconstruct  a  day's  values  of  Sq  for  H,  D,  or  Z 
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The  BAST!'  language  computer  program  using  the  above  formulation  and  Tables  2, 
1,  and  4  to  determine  the  daily  values  o‘  Sq  at  a  given  location  and  date  Is 


given  In  Appendices  A  and  B. 


The  computed  Sq  changes  nake  five  assumptions-  The  first  is  that  the 
four  Fourier  coefficients  (24,  12,  ft,  and  f>  hr)  can  adequately  represent  the 

quiet  daily  variation.  This  convenient  feature  occurs  because  the  principal 
Ionospheric  current  generation  depends  upon  the  diurnal  changes  of 
conductivity  and  winds.  Global  vortex  patterns  of  such  current  and  sudden 
changes  In  the  Ionosphere  at  sunrise  and  sunset  are  mainly  responsible  for  the 
harmonics  at  12,  ft,  and  ft  hr  with  steadily  decreasing  amplitudes.  Although 
higher  harmonics  are  present  at  times,  those  usually  are  hurled  In  a  general 
background  field  level  that  Is  more  like  "noise"  in  character  and  contributes 
to  the  error  of  the  analysis. 

The  second  assumption  Is  that  the  principal  variations  In  the  four 
Fourier  components  throughout  the  year  at  an  observatory  are  the  annual  and 
semiannual  changes.  The  second  Fourier  analysts  provided  values  of  the  annual 
mean,  annual  variation,  and  semiannual  variation  of  the  four  harmonic 
coefficients.  Such  variations  are  expected  because  of  annual  and  semiannual 
changes  that  are  observed  in  the  earth's  heating  and  ionization  during  Its 
yearly  path  around  the  sun.  There  Is  also  a  semiannual,  equinoctial 
interaction  of  the  earth's  magnetosphere  with  the  solar  wind.  Other  variation 
harmonics  have  been  found  In  the  foci  locus.  Longer  perl. id  effects,  such  as 
those  associated  with  the  1  1  - yr  solar  cycle  also  occur.  For  the  short  term, 
one -year  analysts  of  this  study  It  Is  assumed  that  such  changes  are  smaller 
than  the  annual  and  semiannual  variations. 

The  third  assumption  Is  that  the  distribution  of  stations  selected  for 

the  study  ITahle  1}  could  he  used  to  represent  the  smoothed  latitude  change  in 

the  above  coefficients.  Of  course  these  were  the  only  suitable  stations  that 

could  be  found  in  this  region  of  the  earth.  Tt  was  good  to  discover  that  the 

station  coefficients  change  only  gradually  with  respect  to  latitude  (cf. 

Fortunately  sufficient  stations  were  available  to  reflect  the 
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special  Sq  regions  at  the  equator,  mid  latitudes,  auroral  zone,  and  polar  cap. 

The  fourth  assumption  is  that  for  this  North  American  study  the  earth 
induction  characteristics  at  the  important  Sq  frequencies  vary  little  between 
stations  of  adjacent  geomagnetic  latitude  location.  Certainly  for  induction, 
the  longer  period  field  changes  relate  to  deeper  conductivity  regions  of  the 
earth;  and  these  deeper  regions  are  expected  to  change  in  a  more  spherically 
symmetrical  way  than  the  surface  materials.  Observatories  at  the  same 
geomagnetic  latitude  but  differing  in  their  mid  continent,  coastal,  or  island 
location  do  present  conditions  under  which  the  local  induced  currents  will 
contribute  in  different  ways  to  the  locally  measured  field.  The  induced 
currents,  therefore,  are  a  major  source  of  error  in  the  analysis. 

The  fifth  assumption  Is  that  the  Fourier  coefficients  can  be  limited  to 
5°  geomagnetic  latitude  Increments  and  that  the  values  for  any  selected 
observatory's  latitude  may  be  linearly  interpolated  from  the  tables  of  these 
values  (Tables  2,  3  and  A).  This  is  equivalent  to  an  assumption  that  the 
slope  of  the  coefficient  versus  latitude  curve  changes  little  over  a  5° 
Increment.  Figure  13  and  similar  plots  give  some  measure  of  confidence  in 
this  assumption. 

IV.  Results 

The  dally  variations  of  Sq  in  H,  D,  and  7.  for  3°  latitude  increments  on 
representative  equinoctial  and  solstitial  days  (Figure  14)  were  computed  using 
the  full  40  coefficents  of  each  latitude  column  of  Tables  2,  3  and  4..  A  zero 


Figure  14.  Near  here 

level  represents  the  mean  of  the  dtsplaved  variation.  For  this  presentation 
the  distance  between  magnetograms  scales  20  gammas.  Several  general  features 
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should  be  pointed  out.  In  the  display  for  Sq(H)  the  equatorial  region  shows 
the  effect  of  the  very  large  eastward  electroject  current  causing  the  H 
component  to  have  a  range  of  over  1(10  gammas.  There  are  clear  semiannual 
enhancements  in  H  during  the  equinoctial  months  of  March  and  September.  At 
30°  to  40°  we  see  the  lowest  values  of  Sq(H):  this  is  near  the  region  where 
the  ionospheric  Sq  current  vortex  would  be  expected.  At  the  high  latitudes  we 
observe  the  auroral  zone  and  polar  cap  enhancements  with  maximum  values  in  the 
June  solstice  and  continued  large  values  of  field  change  throughout  the  day. 

In  the  Figure  14  display  of  Sq(T>)  note  the  absence  of  extreme 
enhancements  at  the  equator  and  polar  cap  that  were  found  for  Sq(H).  Also 
note  the  regularity  of  the  morning  maximum  and  afternoon  minimum  at  most 
latitudes;  the  transition  (axis  crossing)  shifts  from  near  noon  to  slightly 
earlier  hours  with  increasing  latitude. 

In  the  Figure  14  display  of  Sq(Z)  the  regularity  of  minimum  values  near 
midday  prevails  at  most  mid  and  low  latitudes.  A  clear  semiannual  variation 
occurs  at  the  low  latitudes.  The  field  values  of  Sq(Z)  are  extremely  small 
for  the  polar  cap  winter  solstice  month. 

The  annual  values  of  Sq  are  shown  in  Figure  15.  Recall  that  the  plots 


Figure  15.  Near  here 

are  obtained  from  the  constant  term  of  the  Fourier  series  analysis  for  the 
annual  and  semiannual  variations;  that  is,  rows  I,  2,  II,  12,  21,  22,  31,  and 
32  of  the  Tables  2,  3,  and  4.  This  figure  represents  the  dally  field  levels 
about  which  the  annual  and  semiannual  changes  occur.  It  Is  from  these  plots 
that  the  general  latitude  features  of  Sq  are  described. 

The  annual  and  semiannual  amplitude  variations  of  field  are  represented 
In  Tables  2,  3,  and  4  by  the  coefficients  In  rows  3,  4,  7,  and  8  for  the  12-hr 
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component;  by  rows  23,  24,  27,  and  28  for  the  8-hr  component;  and  by  rows  33, 
34,  37,  and  38  for  the  6-hr  component.  A  plot  of  these  features  at 
representative  latitudes  Is  shown  In  Figure  16.  For  Sq(H)  note  the  well 

Figure  16.  Near  here 

defined  semiannual  change  In  the  24-hr  coefficient  and  Its  harmonics  at  the 
equator.  The  seasonal  changes  of  H  about  disappear  In  the  mid  latitudes  then 
show  a  dominant  annual  change  at  high  latitudes.  At  the  polar  cap  note  how 
the  24-hr  component  dominates  the  amplitude  changes  with  an  annual 
variation.  The  D  annual  change  Is  noticeable  at  all  latitudes.  Here  again 
the  24-hr  component  at  the  polar  cap  dominates  the  picture.  The  significant 
feature  of  Z  Is  the  large  annual  variation  of  the  24-hr  component  at  75° 
latitude.  For  the  more  detailed  latitude  picture  of  the  annual  and  semiannual 
variations  I  have  presented  the  24-  and  12-hr  components  In  Figures  17  and  18. 

Figures  17  and  18.  Near  here 

To  evaluate  the  accuracy  of  the  analytical  representation  of  the  Sq  field 
with  the  limited  number  of  coefficients  I  compared  two  sets  of  values.  The 

first  set,  OBS  (for  "observations"),  was  composed  of  the  hourly  field  values 
scaled  from  the  monthly  quiet-day  curves  representing  the  Sq  field  change  at 
each  station.  As  the  reader  may  recall  these  values  were  the  average  of  the 
month's  quietest  days  with  the  lunar  variations  removed.  To  limit  the  size  of 
the  data  set  only  the  equinoctial  and  solstitial  months  were  used;  this 
selection  gave  96  data  points  for  each  observatory.  The  second  set,  FIT  (for 

"Fourier  fit"),  was  obtained  from  the  analytical  representation  of  the  quiet 

day  field  at  the  same  hours  on  the  16th  day  of  the  selected  months.  All 


observatory  locations  of  Table  1  were  investigated. 

Tigure  19c  illustrates  the  correlation  between  the  two  sets.  A  high 

Figure  19.  Near  here 

correlation  value  means  that  the  shape  of  the  daily  and  seasonal  changes  in 
the  observations  and  the  Fourier  model  match  quite  well;  that  is,  the  relative 
contributions  of  the  Fourier  components  is  representative  of  the  station 
data.  The  correlations  for  the  H  component  were  excellent  except  near  the  Sq 
source-current  vortex  region  between  30°  and  40°  and  near  the  auroral  zone 
between  65°  to  75°.  D-component  values  were  low  near  the  equator.  Z  values 
were  low  at  the  auroral  zone  and  polar  region.  The  three  stations  with  the 
best  correlations  for  all  3  field  directions  were  Paramaribo,  Fredericksburg, 
and  Victoria.  Relatively  good  correlations  were  also  obtained  at  Huancayo, 
Honolulu,  and  Baker  Lake.  The  mean  and  median  values  of  the  correlation 
coefficient,  all  stations  considered,  were  (respectively)  0.907  and  0.936  for 
H;  0.918  and  0.938  for  T);  and  0.809  and  0.959  for  Z. 

To  see  how  the  field  magnitudes  compared  I  will  use  the  term  "error"  to 
represent  the  absolute  value  of  the  residual  difference  between  the 
corresponding  elements  of  the  sets  OBS  and  FIT.  The  mean  value  (MEA)  of  this 
error  for  each  station  is  plotted  In  Figure  19a.  The  standard  deviations,  o, 
of  these  means  were  linearly  related  to  the  corresponding  means  (coefficient 
of  regression  r^-0.93)  such  that  o=0.75  MEA.  The  median  values  of  (MED)  of 
the  errors  were  linearly  related  to  the  mean  (r^=0.97)  such  that  MED=0.86  MEA. 

Figure  9a  shows  that  from  the  equator  to  60°  latitude  the  errors  are 

about  2  gammas  or  less.  At  the  auroral  zone  the  errors  in  H  reach  about  12 

gammas,  the  errors  in  D  reach  about  10  gammas,  and  the  errors  in  Z  reach  5.5 

gammas.  The  Z  errors  are  also  high  in  the  polar  cap.  A  better  representation 

18 


of  chese  station  errors  is  the  percent  error  obtained  from  the  ratio  of  the 
mean  error  to  the  range  of  observed  field  values  at  an  observatory  (Figure 
19b) .  The  largest  percentage  errors  occur  for  Z  between  60°  and  70° 
latitude.  Paramaribo,  Fredericksburg,  and  Victoria  have  the  lowest,  all-over 
percentage  errors  for  all  three  components  combined  whereas,  Huancayo, 
Tatuoca,  Tucson,  and  Baker  Lake  show  reasonably  low  all  all-over  values  of 
less  than  about  10  percent.  The  average  percentage  error  for  all  stations  was 
8.0*  in  H,  5.9*  in  D  and  9.6%  in  Z. 

If  the  values  in  Figures  19a  and  19c  are  considered  the  necessary 
measures  of  the  analytical  representation  of  Sq  then  it  is  possible  to  order 
the  observatory  records  such  that  Victoria  is  best,  Paramaribo  and 
Fredericksburg  are  very  good,  and  Huancayo  and  Baker  Lake  are  good.  For 
future  studies  on  induced  effects  of  Sq  these  five  stations  would  best  be 
used  in  the  analysis. 

If  it  is  assumed  that  the  error  is  partitioned  among  the  various  Fourier 
components  of  the  analytical  representation  in  direct  proportion  to  the 
amplitude  of  the  various  components,  then  Figure  20  shows  how  the  mean  error 

Figure  20.  Near  here 

* 

at  each  observatory  (Figure  19a)  is  divided  between  the  Sq  annual  mean,  annual 

variation  and  semiannual  variation.  Of  course,  such  representation  is  a  bit 

artificial  for  two  principal  reasons.  First,  an  obvious  amount  of  the  error 

comes  from  the  presence  in  the  observations  of  field  changes  with  periods 

other  than  the  four  Fourier  components.  Second,  the  causes  of  annual  and 

semiannual  changes  may  be  different  enough  to  load  the  error  into  one  or  the 

other  of  these  long  period  changes  in  difference  to  their  relative  Fourier 

amplitudes;  such  would  be  the  case  if  solar  activity  affected  the  semiannual 
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change  more  than  the  annual  change-  Intuitively,  It  seems  that  the  longer 
period  changes  are  less  susceptible  to  error.  The  partitioned  errors  give  an 
estimate  of  the  reliability  of  the  seasonal  changes  that  I  personally  think 
over-rates  the  magnitudes.  Figure  20  indicates  that  the  annual  and  semiannual 
variations  have  errors  less  than  about  0.5  gammas  in  H,  D,  and  Z  at  latitude 
below  60°.  At  the  auroral  zone  and  polar  cap  latitudes  the  errors  are  larger 
but,  in  proportion  to  the  range  of  annual  and  semiannual  changes  there,  the 
error  is  about  the  same  as  those  at  lower  latitudes. 

V.  Concluding  Remarks 

The  analytic  representation  of  Sq  for  a  year  of  very  low  solar  activity 
was  obtained  from  monthly  mean  values  of  field  for  selected  quiet  magnetic 
days.  Using  the  diurnal  and  three  higher  harmonic  terms  of  a  Fourier  analysis 
the  coefficients  were  obtained  to  represent  each  Sq  for  the  three  vector 
directions,  for  each  month,  and  at  15  selected  locations  principally  in  the 
North  American  hemisphere.  The  yearly  behavior  of  these  coefficients  were 
then  represented  by  an  annual  mean  plus  an  annual  and  semiannual  variation 
obtained  from  a  second  Fourier  analysis.  Value's  of  all  the  significant  terras 
in  the  analyses  were  then  plotted  as  a  function  of  geomagnetic  latitude  and  a 
smoothing  technique  was  applied  to  obtain  tables  of  40  representative  values 
of  Sq  behavior  for  each  5°  latitude  increment  for  each  fixed  component.  With 
Equation  1  these  tables  were  used  to  reconstruct  a  picture  of  the  daily  Sq 
variations  and  their  seasonal  changes  (figure  14).  Errors  in  the 
representation  were  encountered  primarily  because  of  use  of  a  limited  number 
of  Fourier  coefficients  in  the  analysis,  some  random  variability  in  the  quiet 
natural  field,  smoothing  of  values  from  an  uneven  distribution  of 
observatories,  and  differences  in  the  contribution  of  Induced  currents  to  the 
measured  field  at  each  site.  A  measure  of  these  errors  was  taken  to  be  the 
residuals  between  the  analytical  representations  and  the  observations.  There 


was  a  variation  in  the  error  levels  at  the  various  sites  but  on  average  the 
error  was  less  than  10%  of  the  Sq  amplitude  range  for  the  stations.  There  was 
a  high  value  (0.9)  of  the  correlation  coefficients  between  the  analytic 
representations  and  the  observations. 

A  comprehensive  representation  of  Sq  has  been  presented  earlier  by 
Matsushita  (1967).  His  analysis  was  for  1958,  a  year  of  high  solar- 
terrestrial  activity.  He  divided  the  study  into  three  longitude  zones,  one  of 
which  was  the  American  region.  Although  both  north  and  south  latitude 

locations  were  Included  in  10°  increments,  the  results  above  60°  latitude  had 
been  removed  from  the  presentation.  Seasonal  changes  were  represented  by  4- 
month  group-averages  for  Equinox,  June  Solstice,  and  December  Solstice  (E,  J 
and  D)  months.  In  that  study  the  declination  values  of  D  were  represented  In 
degrees. 

There  are  a  number  of  similarities  and  a  few  differences  in  the  results 

of  the  present  and  the  Matsushita  study.  As  might  be  expected  because  of  the 

difference  In  solar  activity  levels,  the  Matsushita  values  of  Sq  were 

significantly  larger.  The  approximate  source-current  focus  latitude  for  1965, 
represented  by  a  noticeable  transition  in  the  Sq(H)  form,  seems  to  be  near  35° 
latitude  In  summer  and  equinox  but  closer  to  25°  or  30°  in  winter;  this  is 
about  the  same  latitude  as  the  1958  results.  The  shape  (i.e.,  the  relative 
contribution  of  the  Fourier  components)  of  the  Sq  changes  in  the  two  studies 
are  quite  similar  as  are  the  directions  of  change  in  amplitude  between 
seasons . 

The  present  study,  In  addition  to  Its  representation  of  a  very  quiet 
solar  year,  has  too  Important  features.  First,  the  presentation  of  the 
results  above  60°  latitude  shows  the  transition  of  Sq  to  the  large  hlgh- 
latltude  amplitudes  with  major  annual  changes  that  are  largely  representative 
of  the  summer-winter  heating  and  Ionization  differences.  The  significantly 


quiet  solar-actlvlty  year  showed  the  Sq(H)  and  Sq(Z)  transition  to  auroral 
latitude  behavior  at  about  70°.  The  consistency  of  the  Sq(D)  variation 
pattern  at  lower  latitudes  seems  to  be  preserved  at  the  high  latitudes. 

Second,  a  high  resolution  of  the  annual  and  semiannual  changes  Is 
possible  In  the  present  method  of  analysis.  These  changes  were  separately 
plotted  In  Figures  16,  17,  and  18.  Olson  (1970)  had  studied  the  seasonal 
magnetospher 1c  configurations  that  could  contribute  to  a  surface  Sq-type 
variation  and  plotted  these  surface  effects  at  0°,  30°,  and  60°  latitude. 

Comparing  ray  Figure  16  with  Olson's  (1970)  Figures  9  and  10  I  find  that  at  the 
equator  the  form  of  the  two  representations  is  similar  for  H,  only  the  present 
values  are  about  20  times  those  prescribed  by  Olson-  For  equatorial  Z,  Olson 
shows  no  semiannual  terra  and  my  present  annual  value  seems  to  be  in  the  same 
phase  but  about  twice  the  amplitude  of  the  magnetospherlc  model.  The  seasonal 
change  in  the  equatorial  midnight  field  levels  (Figure  II)  follow  the 

magnetospherlc  expectations  (Campbell,  1981)  and  also  do  not  resemble  the 
observed  Sq.  At  30°  latitude  the  presently  reported  annual  variation  in  Sq(H) 
Is  negligible  In  size  whereas  Olson  shows  a  4  gamma  annual  variation  with  a 
summer  maximum.  My  values  for  Sq(Z)  are  larger  and  out  of  phase  with  his  at 
this  latitude.  At  60°  the  phases  of  the  compared  H  variations  are  similar  but 
my  amplitudes  are  about  6  times  those  in  Olson's  figures;  my  Z  component 

values  are  out  of  phase  with  Olson's  model  and  the  amplitudes  also  differ 
considerably.  The  annual  and  semiannual  changes  anticipated  by  the 
magnetospherlc  model  do  not  match  the  observed  changes  sufficiently  to 

encourage  an  assumption  that  the  magnetosphere  Is  more  than  a  minor 

contributor  to  the  Sq  variation. 

Two  final  cautious  reminders  are  appropriate.  First,  the  field  changes 
presented  here  are  determined  for  the  quiet  days  of  a  year  that  was  Itself 
very  quiet  with  respect  to  solar  activity.  The  Ionospheric  conductivity,  wind 
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and  tidal  systems,  and  location  of  the  auroral  zone  should  be  expected  to  be 
different  on  days  of  higher  solar  activity.  Second,  the  separation  of  the 
induced  Internal  field  from  the  external  source  contributions  has  not  been 
Included  In  the  present  study;  such  separation  Is  planned  for  the  future. 

The  present  results  represent  the  best  determination  of  the  quiet 
geomagnetic  variations,  Sq,  within  North  America  to  be  expected  as  minimum 
field  variation  levels  under  very  quiet  conditions.  The  results  Illustrate 
the  suitability  of  four  Fourier  components  to  describe  the  Sq  and  to  detail 
the  annual  and  semiannual  changes  of  these  components. 
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Appendix  A 


BASIC  Program  for  storage  of  Tables  2,3,  and  4  on  tape  file. 


1  DISP  "PRuGRAM  ' OAFILE'  TO  |_n 
HO  FILE  ONTO  TAPE" 

2  REM  * F 0 R  USE  IN  CREATING  FIL 
E  0 F  S Q  F 0 U R I E R  C OEFFICIENT S 

US  LATITUDE:*: 

3  DISP  "INPUT  NEW  FIlE  ARRAY  N 
AME  AS  FT" 

4  INPUT  FT 

5  DISP  "INPUT  MUMPER  OF  ARRAY 
COLUMNS  C” 

6  INPUT  C 

7  DISK  “INPUT  NUMBER  OF  ARRAY 
ROWS  R “ 

8  INPUT  R 

9  DISP  “BEFORE  RUN  FIX  90  DIM 
W  <  R  >  “ 

30  CREATE  F*,C,3*R 
40  REM  t  F  I L  E  C  C  0  L  <  R  E  C  O  p  D  S  .  .  R  R 
OWSC ITEMS  EA >  WITH  8  BYTES  P 
ER  ROWf. 

50  OPTION  BASE  1 
60  REN  * A BODE  TO  START  COUNT  AT 
1  RATHER  THAN  0:*: 

70  ASSIGN#  1  TO  FT 

80  REM  *  ABQOE  FOR  BUFFER* 

90  DIM  W < 4 0 > 

100  FOR  L=i  TO  C 
110  DISF . 

120  DISP  “NEXT  R  VALUES  FOR  COLU 
MN  “  ;  L 

130  FOR  1=1  TO  R 
140  INPUT  W(I) 

150  NEXT  I 

160  BEEP 

161  REM  *  BtEP  FOR  END  OF  DATA  E 
NT  R 

1  7  0  0  I  S  P  “  V  A  L  U  E  S  F  0  R  C  0  L  U  M  N  “  L 
1  §  0  F  0  R  I  =  i  T  Cl  R 

190  DISP  “H“  I  ;  “  =  “  ;  1-K  I  ) 

191  REM  *  CORRECT  W3  =  2  WITH  W ■,  3  '•  = 
3  ON  0  I  SPLAY ; PRESS ' END  LINE' 

* 

200  NEXT  I 

210  DISP  "CHECK, CORRECT,  THEN  PR 
ESS  COnT” 

220  PAUSE 

230  PRINT#  1  ,  L  *  WO 
2 4 0  P EM  YAP 0 E  LISTS  ALL  PEC 0 R 0 S 
ON  TAPE* 

250  DISP  “END  R  ENTRIES" 

260  NEXT  L 

270  ASSIGN#  1  TO  * 

230  REM  * ABOVE  CLOSES  BUFFER* 

2 P  0  D 1 8 P  " E  N  U  FILE  X  D A T  0  A P  k A  V " 
3O0  END 


Appendix  B 


BASIC  program  for  plotting  of  Sq  from  tape  file  of  Tables  2,3,  and  4. 


IS  REM  *  SOP  LOT*' 

20  DEF  FNR <  02  >  =  I  NT < D2* 1 90+ . 5 > 
✓  1  00 

30  REM  * ABOVE  POUND IMG  FCN  FOR 
2  PLACES* 

40  DISP  "ENTER  STATION  NAME" 

50  INRuT  Nf 

51  PRINT  Nf 
60  DEG 

70  DISP  "ENTER  GEOGRAPHIC  LAT  1 
N  DECIMAL  DEG  WITH  NORTH  ♦ 

80  INPUT  G 1 

90  PRINT  Gli-OEG  LATITUDE" 

100  DISP  “ENTER  GEOGRAPHIC  LONG 
IN  DECIMAL  DEG  WITH  EAST  +  " 
110  INPUT  C2 

120  PRINT  G2  ’ UEG  LONGITUDE" 

.130  C1=C0SC  1  1 . 3> 

140  C  2  =  S I  N  C 1 1  3  > 

150  C3=C0S  v  G 2 - 2 8 9  .  6 > 

1  60  C 4  =  S  I N  ■:  G2-289  6  > 

170  L1=ASNCSINCG1  > tC  1  +  C0S< G 1  >*C2 
*C3> 

,180  L2  =  ASN  <  COS  <  G 1  >  *C4-'C0S<  LI  >> 

181  IF  L £ >  =  0  THEN  104 

182  L2--360  +  L2 
1 J3  GOTO  190 

184  IF  L 2  <' 3 6 0  THEN  ISO 

185  L 2  =  L 2 - 3 6 O 

190  PRINT  FNR':  LI  ■ ;  “DEG  GEOMAGNET 
IC  LATITUDE" 

200  DISP  FNR CL  1 > ; "DEG  GEOMAGNET I 
C  LATITUDE" 

210  PRINT  FNRCL2>; "DEG  GEOMAGNET 
IC  LONG" 

220  DISP  FNRCL2); "DEG  GEOMAGNET  I 
C  LONG" 

221  IF  L 1 <0  OR  LI >90  THEN  223 

222  GOTO  230 

223  BEEP 

224  DISP  "GEOMAG  LAT  OUTSIDE  RAN 
GE  -SELECT  NEW  LOCATION  AND 
RUN" 

225  PRINT  "NO  COMPUTATION" 

226  PRINT  "  " 

227  PAUSE 

230  DISP  "SELECTED  ANALYSIS  DATE 

It 

231  DISP  " DAY= " > 

232  INPUT  D 

233  DISP  "MONTH*"; 

234  INPUT  M 

237  Y= 1 965 

238  IF  M= 1  OR  M=2  THEN  242 

239  Y 1 =  Y 
249  M 1 =  M  +  1 

24 1  GOT  0  244 

242  Y 1 =  V - 1 

243  M 1 =  M ♦ 1 3 

244  1 1  =D- I  NT  c  .  75* I  NT  (.  Y  1  /  1 00  >  -7  >  + 
f  NT  <.  365  254  Y 1  >  +  INTC30 . 6*M1  > 


245  J2=  1  -  I N  T  v  75*  I  NT  (  •.  Y-l  >,'100  >- 
7>  +  I  NT  (.365  25*  t  V  -  1  )  +  I N  T  c  30 
6*14' 

246  D 1  =  J 1 - J2+ 1 

250  PRINT  “DAY=" . D . "MONTH*" . M, " D 
AY  NUMBER  FOR  YEAR  =  ",D1 
260  M  1  =  '  0  1  -  1  >  3 Pi  4  16  7 
280  DISP  FNR UU  /  "  IS  ANALYSIS  MO 
NTH'* 

290  U 1  =24-G2-'  1  5 

291  IF  U 1 <24  THEN  300 

292  U 1  =U  1 -24 

300  PRINT  "UT  OF  LOCAL  MDT  =";FN 
R(U1  > 

310  DiSP  FNPt U1 > IS  UT  AT  LOCAL 
NOON " 


320 

REM  *N=t, 
4  ,  12-8 ■ 6 

HP 

3.4  FOR  PERIOD 
WHERE  T=24 'N" 

-■  2 

330 

REM  *P=1, 
SE  •:  2  '  * 

2' 

FOR  AMP-:  1  )  OR 

PHA 

348 

350 

REM  *0=1, 
I  ANN  < 2  ■ 

L  3  =  L  1  s  5  M 

* 

FOR  ANN'-  1  )  OR 

SEM 

360  REM  *  SELECT  BRACKETING  L* 

370  FOR  L  = 1  10  19 

380  IF  L3>=l  AND  L-3CL+1  THEN  400 
390  NEXT  L 

400  DISP  "LOWER  TABLE  LATITUDE"- 
5*:  C  L  -  1  >  ,  "DEG" 

410  OPTION  BASE  1 

411  FOR  A- 1  TO  3 

412  IF  A= 1  THEN  430 

413  IF  A=2  THEN  432 

414  IF  A=?  THEN  434 

430  ASSIGN#  1  TO  "SOHLAT" 

431  C  f  =  "  H  "  G 0 TO  44 0 

432  ASSIGN#  1  TO  "SODLAT" 

433  Cf="0"  i?  GOTO  440 

434  ASSIGN#  1  TO  "SOZLAT" 

435  Cf=“Z"  A  GOTO  440 
44  0  DIM  W •:  4 O  ’ 

441  DIM  H 7 1  4 0 '> 

442  DIM  Y 1  4 0  > 

443  DIM  HC19.401 
450  FOR  X=L  TO  L+l 

460  READ#  1  ,  X  ;  W (  ) 

461  DISP  "VALUES  AT  LAT "  ;  FNR •' 5* v 
X-i>> 

470  FOR  1=1  TO  40 
480  H  <  X  .11=  W  C  I  > 

490  DISP  HCX-1) 

500  NEXT  I 
510  NEXT  X 

550  FOP  X=1  TO  40 

551  y <. x  > = h •: l + 1  ,x  '-h-:l,x> 

560  IF  X=2  OR  X= 1 2  OP  X=22  OR  X= 
32  THEN  580 
570  GOTO  640 

580  IF  Y <  X  >'  =  12  THEN  61 0 
590  IF  Y < X >  <=— 1 2  THEN  630 


60O  GOTO  *31 
610  Yl;0*YfX  >-24 
620  GOTO  *31 

630  Yv X >= Y<X >+24 

631  H3'  X  >  =  ■'  L3-L  >  *  V  O' >  +  H «  L ,  X) 

632  IF  H3  <  X  > ' 0  THEN  830 

633  Hi' X  *=H3' X  >  +  24 

634  GOTO  830 

640  IF  X-4  OR  X  =  6  OR  X  = 1 4  OR  X=  1 
6  OR  X  =  24  OR  X=2*  0R  X  =  34  OF: 
X  =  36  THEN  6* O 
6*50  GOTO  730 
660  IF  Y  YX>:  =6  THEN  *80 
670  IF  Y'X>Y=-6  THEN  710 
680  GOTO  711 
680  Y  Y  X  1  =  Y  Y  X  >  -  1 2 
700  GOTO  711 
710  Y  Y  X  >  =  Y  Y  X  '•  +  1  2 
71  1  H  3  Y  X  >  =  Y  L  3  -  L  >  *  Y  ■  X  1  +  H  >  L  ■  X  > 

712  IF  H  3  Y  X  >  >  0  THEN  830 

713  H  3  <  X  >  =  H  3  Y X >  +  12 

714  GOTO  870 
720  GOTO  820 

730  IF  X=8  OR  X =10  OR  X=1S  OR  X= 
20  OR  X  =  2 8  OR  X  =  30  OR  X=38  0 
R  X=40  THEN  750 
740  GOTO  S20 
750  IF  YCX)>«3  THEN  780 
760  IF  Y  Y  X  )  <  =  -  3  THEN  800 
770  GOTO  801 
780  Y<X>=Y'X>-6 
790  GOTO  801 

800  Y  Y  X  >  =  Y '  X  >  +  6 

801  H  3  Y  X  >  =  Y  L3-L  >  *Y  Y  X  >  +  H  Y  L  .■  X> 

802  IF  H3 < X > > 0  THEN  830 

803  H 3  (  X  >  =  H 3  <  >0+6  • 

804  GOTO  830 
810  GOTO  820 

820  H  3  Y  X  >  =  Y  L  3  -  L  >  *  Y Y X  >  +  H  *  L  , X) 

830  REM  *ABOVE  GIVES  I NTERP0L8TE 
D  VALUES* 

840  0 1 SP  "  TAB  VAL"  ;  FNF:  Y  X  FN 

R  Y  H  3  Y  X  :>  :> 

850  NEXT  X 

860  ASSIGN#  1  TO  * 

870  REM  *NOU  HAVE  40  INTERPOLATE 
0  VALUES  OF  TABLE  FOR  GEOMAG 
LAT  LI# 

880  REM  ♦  N = P  E  R I O 0  SELECT  SO  T  =  24 
✓N* 

831  REM  *R=1  FOR  AMP;  2  FOR  PHA3 
E  * 

:  82  REM  *0=1  FOR  ANN;  2  FOR  SEMI 
ANN* 

883  REM  * M 1 =  MONTH  VALUE* 

384  REM  *F - N , R  >  AND  F 1 <  0 > 

920  DEG 

930  DIM  FY4,2> 

940  DIM  F 1 Y 2 ) 

950  FOR  N= 1  TO  4 


960  FOR  F:  =  1  TO  2 
970  F OF  0=1  TO  £ 

980  Rl=2*R-2 

990  F 1 Y  Cn  =H3 '  3  +  R 1 +4*  <  0-1  >  + 1  0 *  v  N- 
1  >  >  *  S  I  N 1  30*0* < M J -H3  <  4  +  R 1 +4  * Y 
Q- 1 >  +  10*<  N-l ) ) >  +  90; 

1000  NEXT  0 

1010  FyN • P>=F1 <  1 >+Fl Y2> 

1011  DISP  "  F  y  “  ;  N  ;  R ;">=“; F NR  f 

F  Y  N  •  R  >  > 

1020  NEXT  P 
1030  NEXT  N 
1040  DIM  H  4  Y  2  4 1  ) 

1050  DIM  Y  <■  4  :• 

10*0  FOR  N= 1  TO  4 
107O  FOR  R= 1  TO  2 

1080  F  (  N  •  F‘  >  =  F  v  N  ,  F,‘  >  +  H  3  Y  R  +  1  0  :*  >;  N  -  1  > 

> 

1090  DISP  "TOT  F< " ; N; ", " ; R; “ >  =  " , 
F  NF'  v  F  N  .  R  >  > 
lioo  next  r 
1110  next  n 

1120  FOR  T 1 = 1  TO  £41  STEP  1 

1121  T  =  ' T 1  - 1  ■ ' \ 0 
1140  FOP  N= 1  TO  4 

1 1 50  M  H ' =F ' N •  1 > tS I N  Y 1 5*  N*  Y  T - F  Y  N 
, £ }  )  +90  > 

1160  NEXT  N 

1170  H4  <  T 1  >  =K '  1 >  +K '  2  >  +  P ''  3  '  +  M  4  ■ 

1171  DI  SP  C *  ;  “  Y  "  ;  T  ;  "  >  =  "  ;  F N R  1  H 4  T 
1 >  > ;  "GAMMAS" 

1180  NEXT  T 1 

1181  BEEF' 

1182  PRINT  "  "  9  PRINT  " "  6  PRINT 

M  II 

1190  FOR  T  = 1  TO  24 

1191  T 1  =  1 0  *  T  + 1 

1  20O  PR  I  NT  C  l  ;  "  <  "  ;  T  ;  "  >  =  "  ;  FHR  •  H4  Y 
T 1 >  > ; "GAMMAS " 

1210  NEXT  T 

1211  PRINT  . .  PRINT  "" 

1212  PRINT  C 1'  .  "  "  ;  N*  .  "  ","DaV 

; D; "Month" ; M 

1213  PRINT  "  AMPLITUDE  •' G 

a mm as >  " 

1214  PRINT  " -20  -10  0 

10  20" 

1220  PEN  1  9  GCLEAR  '  START  GRAF* 
H1CS  SCALE 

1230  V 1 =-20  '  BOTTOM  SCALE  SET 
1240  . 2=20  '  TOP  SCALE  SET 
1250  SCALE  0 .24, Y1  •  V 2 
1260  XAXIS  Y 1,1,0. 24  9  XAX I S  Y2, 
1,0,24 

1270  XAXIS  < Y 1 +  Y2 >^2,0  9  XAXIS  Y 
Y  1  +  Y  2  >  /  2  +  1  0 , 0  9  XAXIS  Y  1  +  V 
2)^2- 10,0 

1280  YAXIS  0 , 2 , Y 1 , Y 2  9  YAXIS  24, 
2, Y1 , V  2 

1290  FOR  X= 1  TO  3 


example  -> 


l  ?i?0 

1310  ! 

1311  f 

1312 
1320  I 
1330  : 

1340 

i 

1341 

1342 

1343 

1344  i 
134  5 
134E:  , 

1350 

1351 

1352 

1354 

1355 

135 6 

1357 
135S 

,  1  3*50 

t361 
1  -62 
1363 
1  7:?0 
13:51 
1  3  9  O 
1  391 
1  396 
1397 

1399 

1400 
1410 
1440 
14  50 

1451 

1452 
1460 
1  170 
1  4  SO 

14  90 
1500 
1510 
1520 
1530 
1540 

1550 
1560 

15  .. 

50 

1551 
1582 
1590 
1600 
1610 
1620 


YRM5  6a: 
HEM  • 
PEMUP 
Z  =  0 

FOP  T 1 = 1 

X=M  1-1  * 

IF  Y<=20 


: . O . y 1 • Y2 

E HO  S C O L E  ‘SET 


TO  241  STEP  1 
■  i  ft  i»  Y  =  H4  •  T  1  1 
RNO  Y .*=-20  THEM 


IF  Y  >  20  THEM  1343 
IF  V  <-20  THEM  1356 
IF  Z=0  THEM  1345 
OO TO  135 0 
PEMIJP 
2-1 

"iF*  Y  >  20  HMD  Y  060  T1 
IF  2=1  THEM  1354 
GOTO  1355 
PEMUP  @  2  =  0 
GOTO  1390 
IF  2=0  THEM  1358 
GOTO  1360 
PEMIJP  0  2=1 
IF  YC-20  HMD  Y>=-60 
81 

IF  2=1  THEM  1363 

GOTO  13:50 

PEMUP  8  2=0 

PEN  1  ii  Y  =  Y  +  :50  i5  GO 

PEN  1  'i  Y  =  Y  +  40  <i  GO 

PEN  1  Y  =  V - 8 0  8  GO 

PEN  1  i?  Y  =  Y - 4 0  it  GO 

IF  2=0  THEM  1399 

PEMUP  it  2  =  0 

PEN  1 

plot  h.y 

NEXT  T 1 

PEMUP  1  STRPT  LRBEL 
LD1P  O  it  FEN  1 
HOVE  .2.18.1 
LRBEL  CS 
MOVE  1.-19  9 
LRBEL  “0“ 

HOVE  6  1,-19. 9 
L  R  B  t  L  "  6  “ 

HOVE  12  1,-1 9. 9 
LRBEL  “12" 
hO'-’E  18  1,-1 9. 9 
LRBEL  "18" 

HOVE  22  8, -19. 9 
LREEL  "24“ 

COPY  8  NEXT  R 
CUSP  "IF  WRNT  NEXT 
ESS  CONT" 

FOP  1=1  TO  50 
PEEP  I tRNO+l , 50 
NEXT  I 
l  PHUSE 

M 1  =  Ii  1  +  1  H  M  =  f1  + 1 
»  GOTO  411 
I  END 


THEN  1391 


-60  THEM  13 


GOTO  1400 
G  0  TO  14  O  u 
GOTO  1400 
GOTO  1400 


NEXT  MONTH  PR 


15  0  E  G  G  E  u  1 1  h  1 , 

4  7  06 i.m  I  I  TuCF 

27u  ['[  ,  I  i  UL'F 

15  DC 6  'jL'.'I’h  ,t(E  V  !  C  LmTIIIJOE 
339  7  4  DE'j  ijE  OHhGM t  T  I  C  uONu 
DRV  =  M  H *.* ;  i  ;  H  =  3 
DRY  MU  HE  L  F  F  *.*  F  YE  HP  =  80 
Li  T  uh  L'.-i  HL  HC'l  =  6 


HC  1  ■  = 

H  ■.  2  .  = 
H  (.  j.  •  = 
H  *.  4 
H<  t  >  = 
H  (  6  >  - 
H  (  7  1  = 
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Figure  Captions 

1.  The  daily  variation  in  declination  Sq(D),  the  solar  zenith  angle  X,  and 
the  Chapman  function  /cos  X  for  Honolulu,  Hawaii  on  quiet  days  in  March 
1965.  The  field  variation  is  displayed  about  its  mean  position  of 
11°34.56'  E.  Local  noon  occurs  at  about  22:32  UT  here. 

2.  The  spectral  composition  of  daily  variations  in  declination  of  Tucson, 
Arizona  for  all  days  of  March  1965.  a)  The  spectra  between  28  days  and  6 
hr  period;  b)  The  spectra  between  6  hr  and  5  min  period. 

3.  The  monthly  average  of  quiet  day  field  values  when  K_<2+  in  1965  for  the  H 
component  of  field  at  Fredericksburg,  USA;  Huancayo,  Peru;  and  Toolangi, 
Australia.  The  scale  to  the  left  should  be  multiplied  by  6  to  give  the 
Huancayo  amplitudes. 

4.  Difference  in  appearance  of  magnetogram  H-component  values  on  quiet  (21 
June  1965,  Ap  index  -  2)  and  active  (15  June  1965,  Ap  *  19)  days.  Station 
letter  codes  are  given  to  left  with  the  corresponding  geomagnetic 
latitudes.  Note  the  change  in  amplitude  scales  for  the  station  values. 

5.  Results  of  a  superposed-epoch  analysis  of  March,  1965,  Honolulu 

declination  variation  residuals  in  lunar  time  as  described  in  the  text. 
Note  the  large  semidiurnal  variation  of  the  lunar  tidal  component  L.  0000 
lunar  time  corresponded  to  0450  UT  on  March  3,  1965.  Zero  on  the 

amplitude  scale  corresponds  to  11°34.56/  East  declination. 

6.  Comparison  of  computed  lunar  field  component  to  the  magnetogram  record  of 
Honolulu  declination  (D)  component  of  field  on  March  20,  1965.  Zero  on 
the  amplitude  scale  corresponds  to  11°34.56'  East. 

7.  Variation  of  Chapman  factor  with  geographic  latitude  and  season.  Dashed 
lines  represent  Tropic  of  Cancer  and  Arctic  Circle. 


8.  Daily  variation  in  Universal  Time  of  solar  zenith  angle  (solid  curves) 
and  Chapman  factor  (dashed  curves)  for  summer  solstice  (SS),  equinox  (E), 
and  winter  solstice  (WS)  at  Boulder,  Colorado. 

9.  Zonal  mean  meridional  circulation  in  the  earth's  thermosphere  during 

equinox  and  solstice  for  three  levels  of  auroral  activity  (Roble,  1981). 

10.  Magnetospheric  contours  of  AB  the  variation  of  the  observed  field  from  a 

main  field  model,  in  the  geomagnetic  noon-midnight  meridian  plane  In 

quiet  conditions  (Suguira  et  al*,  1971). 

11.  Seasonal  variation  of  midnight  levels  of  H,  D,  and  Z  field  components  for 
5°  increments  from  the  equator  to  30°  geomagnetic  latitude  (Campbell, 
1981) . 

12.  Block  diagrraas  of  the  contributors  to  ionosphere  current  and  its  surface 

observation.  a)  general  diagram,  b)  details  of  upper  part  of  general 

diagram,  c)  details  of  lower  part  of  general  diagram. 

13.  Values  for  annual  mean  phase  in  hour  of  first  amplitude  maximum,  for  12- 
hr  component  of  Sq(H).  Point  values  are  indicated  with  the  associated 
station  codes  (Table  1)  and  line  segments  connecting  these  values. 
Smooth  curve  is  obtained  from  5-point  smoothing  of  values  read  from  line 
segments  at  2.5°  intervals  with  natural  extensions  for  the  curve  ends. 

14.  Daily  variation  of  Sq  in  5°  geomagnetic  latitude  Increments  for  March, 
June,  September  and  December.  Scale  to  right  indicates  20  gammas  between 
rows,  a)  Sq(H);  b)  Sq(D);  c)  Sq(Z). 

15.  Annual  mean  values  of  Sq  in  5°  geomagnetic  latitude  increments  for  H,  D, 
and  Z.  Scale  to  right  indicates  25  gammas  between  rows. 


16.  Annual  and  semiannual  amplitude  variations  of  Sq  for  the  separated  24-, 
12-,  8-,  and  6-hr  component  periods.  Representative  plots  at  0°,  30°, 
60°,  75°,  and  90°  geomagnetic  latitude  are  shown  for  the  H,  D,  and  Z  field 
directions. 

17.  Annual  and  semiannual  variation  In  24-hr  component  amplitude  of  Sq  for  5° 

geomagnetic  latitude  Increments.  Scale  to  right  Indicates  5  gammas 

between  rows  In  the  H,  D,  and  Z  field  direction  columns. 

18.  Annual  and  semiannual  variation  In  12-hr  component  amplitude  of  Sq  for  5° 

geomagnetic  latitude  Increments.  Scale  to  right  Indicates  5  gammas 

between  rows  In  the  H,  D,  and  Z  field  direction  columns. 

19.  Measure  of  fit  of  Fourier  analysis  values  (FIT)  to  station  observations 

(OBS)  using  hourly  values  In  March,  June,  September  and  December,  for  15 
stations  plotted  at  their  respective  geomagnetic  latitudes.  a)  (bottom) 
mean  residual  (error)  values  as  difference  between  FIT  and  OBS;  b) 
(middle)  residuals  as  percent  of  range  of  OBS;  and  c)  (top)  correlation 

coefficients  between  FIT  and  OBS. 

20  Estimated  error  between  Fourier  analysis  values  (FIT)  and  station 

observations  (OBS)  partitioned  between  the  annual  mean,  annual  variation, 
and  semiannual  variation  in  proportion  to  the  respective  Fourier  component 
amplitudes.  Values  for  15  stations  plotted  at  their  respective 

geomagnetic  latitudes  Is  illustrated  for  a)  (top)  Sq(H),  b)  (middle) 

Sq (D) ,  and  c)  (bottom)  Sq(Z). 
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